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Abstract 

A set of three-dimensional failure criteria for laminated fiber-reinforced composites, 
denoted LaRC04, is proposed. The criteria are based on physical models for each 
failure mode and take into consideration non-linear matrix shear behaviour. The 
model for matrix compressive failure is based on the Mohr-Coulomb criterion and it 
predicts the fracture angle. Fiber kinking is triggered by an initial fiber misalignment 
angle and by the rotation of the fibers during compressive loading. The plane of 
fiber kinking is predicted by the model. LaRC04 consists of 6 expressions that can 
be used directly for design purposes. Several applications involving a broad range 
of load combinations are presented and compared to experimental data and other 
existing criteria. Predictions using LaRC04 correlate well with the experimental 
data, arguably better than most existing criteria. The good correlation seems to be 
attributable to the physical soundness of the underlying failure models. 


1 Introduction 


The greatest difficulty in the development of an accurate and computationally 
efficient numerical procedure to predict damage has to do with how to model 



the material micro-structural changes and how to relate those changes to the 
material response. The results of the Word Wide Failure Exercise (WWFE) 
[1] indicate that — even for simple unidirectional layup and in-plane loading — 
current failure models and criteria fall short of providing a satisfactory de- 
scription and prediction of failure. In fact, the mechanisms that lead to failure 
in composite materials are not fully understood yet. 

If composite materials are to be used in structural applications, then the un- 
derstanding of how each failure mode takes place — i.e. having a physical model 
for each failure mode — becomes an important point of concern. These phys- 
ical models should establish when failure takes place, and also describe the 
post-failure behaviour. For instance, a physical model for matrix compression 
failure should predict that failure occurs when some stress state is achieved, as 
well as what orientation should the fracture plane have and how much energy 
should the crack formation dissipate. 

The main failure modes of laminated fiber-reinforced composites are: 

Delamination. Composite materials made of different plies stacked together 
tend to delaminate. The bending stiffness of delaminated panels can be signif- 
icantly reduced, even when no visual defect is visible on the surface or the free 
edges. The physics of delamination is to a certain degree understood, and one 
of the best numerical tools to predict the propagation of delamination consists 
on the use of Decohesion Elements. These elements have been developed [2,3], 
and implemented in a commercial Finite Element (FE) code [4]; no further 
discussion of delamination as a failure mode is thus carried here. 

Matrix compression failure. What is commonly referred to as matrix com- 
pression failure is actually shear matrix failure. Indeed, the failure occurs at 
an angle with the loading direction, which is evidence of the shear nature of 
the failure process. 

Fiber compression failure. This failure mode is largely affected by the resin 
shear behaviour and imperfections such as the initial fiber misalignment angle 
and voids. Typically, kinking bands can be observed at a smaller scale, and 
are the result of fiber micro-buckling, matrix shear failure or fiber failure. 

Matrix tensile failure. The fracture surface resulting from this failure mode 
is typically normal to the loading direction. Some fiber splitting at the fracture 
surface can usually be observed. 

Fiber tensile failure. This failure mode is explosive. It releases large amounts 
of energy, and, in structures that cannot redistribute the load, it typically 
causes catastrophic failure. 

In this work, the LaRC03 plane stress criteria [5] are extended to account for 


2 



general three-dimensional (3D) loading and for in-plane shear non-linearity. 
Matrix compressive failure is addressed with a modified 3D version of the 
Puck [6] matrix compression failure criterion. For the fiber failure in compres- 
sion, a 3D criterion based on Argon’s [7] approach considering matrix failure 
prompted by material imperfections is used in a framework similar to the one 
proposed in LaRC02 [8] and LaRC03 [5] . A failure model for matrix in tension 
and shear is derived from Dvorak and Laws’s [9] fracture mechanics analyses 
of cracked plies, as in LaRC03 [5]. 

Section 2 presents the background on damage mode-based failure criteria. The 
LaRC04 failure criteria are developed and proposed for each failure mode in 
Section 3. Examples and validation follow in Section 4, and conclusions are 
presented in Section 5. A summary of the LaRC04 criteria is presented in 
Appendix. 

In the following sections, the index 1 refers to the longitudinal (fiber) direction, 
the index 2 refers to the in-plane transverse direction and the index 3 refers 
to the through-the-thickness direction. 


2 Damage mode-based failure criteria 


Strength-based failure criteria are commonly used to predict failure events in 
composite structures. A large number of continuum-based criteria have been 
derived to relate stresses and experimental measures of material strength to 
the onset of failure. Paris [10] discusses the ad hoc nature of the formulation 
of most strength-based criteria. 

Hashin [11,12] and Puck [13,14] are credited for establishing the need for fail- 
ure criteria that are based on failure mechanisms. In his 1973 proposal [11], 
Hashin used his experimental observations of failure of tensile specimens to 
propose two different failure criteria, one related to fiber failure and the other 
related to matrix failure. The criteria assume a quadratic interaction between 
the tractions acting on the plane of failure. In 1980 [12], he introduced fiber 
and matrix failure criteria that distinguish between tension and compression 
failure. Given the difficulty in obtaining the plane of fracture for the ma- 
trix compression mode, Hashin used a quadratic interaction between stress 
invariants. Such derivation was based on logical reasoning rather than mi- 
cromechanics. Although the Hashin criteria were developed for unidirectional 
laminates, they have also been applied successfully to progressive failure anal- 
yses of laminates by using in-situ unidirectional strengths to account for the 
constraining interactions between the plies [15] . The two-dimensional versions 
of the failure criteria proposed by Hashin in 1973 and 1980 are summarized 
in Table 1. 
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Table 1 

Hashin criteria for plane stress 


Constituent Tension 


Matrix F/ M = (^) 2 +(ff) 2 

Fiber FI f = (^ p ) 2 + 


FI m 

FI m 


Compression 

= m 2 +m 2 

/ /To \ 2 


Numerous studies conducted over the past decade indicate that the stress 
interactions proposed by Hashin do not always fit the experimental results, 
especially in the case of matrix or fiber compression. It is well known, for in- 
stance, that moderate transverse compression (022 < 0) increases the apparent 
shear strength of a ply, which is not predicted by Hashin’s criterion. In addi- 
tion, Hashin’s fiber compression criterion does not account for the effects of 
in-plane shear, which significantly reduces the effective compressive strength 
of a ply. Several researchers have proposed modifications to Hashin’s criteria 
to improve their predictive capabilities. 

More recently, the WWFE [16] was conducted to assess the real predicting 
capability of the current available failure criteria. Leading researchers in fail- 
ure of composites were invited to participate in a round-robin in which they 
presented their approaches and predictions. 

In the exercise, Hart-Smith [17,18] presented the original version of the max- 
imum strain criterion, as well as a truncated form and a generalized form of 
the criterion. Gotsis et al. [19] used the maximum stress criterion, superposed 
with a modified distortion energy (quadratic polynomial) criterion. McCart- 
ney [20] applied the principles of mechanics at the microstructural level to 
predict damage formation. Rotem [21] used a criterion originally published in 
1973 [11], but with the matrix failure criterion modified in order to account 
for axial stresses. Surprisingly, Sun and Tao [22] used the Hashin-Rotem crite- 
rion [11], even though Sun had proposed previously [23] a criterion for matrix 
cracking that is acknowledged to represent better the matrix failure [8]. Liu 
and Tsai [24] used the Tsai-Wu [25] failure criterion. Wolfe and Butalia [26] 
used a strain-energy based failure criterion, containing a sum on exponents of 
the longitudinal, transverse and shear strain energies. Edge [27] used a phe- 
nomenological approach based on the stress interaction within each failure 
mode, with some similarities to the maximum stress, the Hashin-Rotem [11] 
and Hashin [12] criteria. Zinoviev [28] used the simple maximum stress crite- 
rion. Puck and Schiirmann [6] were perhaps the authors who achieved better 
agreement between their predictions and the experiments. Their criterion is 
phenomenological, as different failure modes are considered. Arguably, their 
model for matrix compression failure possesses the most sound physical basis 
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of the theories proposed in the exercise. 


Several lessons can be learned from the WWFE. Firstly, most criteria were 
unable to capture some of the trends in the failure envelopes of the experimen- 
tal results. Secondly, on what concerns phenomenological failure criteria, most 
expressions proposed to predict each failure mode are still to some extent em- 
pirical. It is somewhat difficult to choose between the criteria due to the lack of 
experimental data needed to validate them against each other. Despite several 
efforts to develop sound phenomenological criteria, non-phenomenological cri- 
teria like Tsai-Wu [25] are often better prediction tools than some phenomeno- 
logical criteria [24], Although test results are not provided in the WWFE 
for several stress combinations that remain open for discussion, significant 
progress was made. From the limited predictive capabilities of the most accu- 
rate analyses available, it is clear that further developments in failure model 
theories and criteria are required before any analysis approach can be used 
with confidence to predict the strength of a typical aerospace composite com- 
ponent. 


3 LaRC04 criteria for matrix failure 

3. 1 Tensile matrix failure 

3.1.1 Critical energy release rates for non-linear shear behaviour 

This subsection generalizes the expression given by Laws [29] for the energy 
release rate for an elliptical crack in a composite with a non-linear shear be- 
haviour. The procedure presented here uses Eshelby’s [30] application of the 
eigenstrain problem to solve the stress field around an oval crack, in a frame- 
work similar to Laws’s [29]. Since the referred procedure is considerably com- 
plex, and yet not devoid of approximations, an alternative generalization that 
leads to the same result is presented in the Appendix. This alternative gen- 
eralization is much simpler, though cruder; it is pragmatic in the sense that 
no complex mathematical manipulations are required, and the result follows 
from the concept of critical energy release rate and intuitive approximations. 


The Eigenstrain problem Consider an infinite solid that undergoes some 
physical process that will generate some inelastic strain. The physical process 
could for instance be plastic deformation or phase transformation. Probably 
due to the later, this inelastic strain is commonly called transformation strain. 
The problem consists in determining the stress and strain fields due to the 
transformation strains. 
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The total strain tensor £ in the transformed region is the sum of an elastic 
strain tensor £ el with the transformation strain tensor £ t \ 


£ = £ el + £ t . ( 1 ) 
The stresses result from the elastic component of the strain, 

a = C-.£ el = C:(£-£ t ). (2) 

Neglecting body forces acting on the body, the equilibrium equation is 

V • a = C : (V • £ - V • £*) = 0 (3) 

and can be written as 

C : V • £ — b* = 0 (4) 

where b* is a fictitious body force distribution defined as 

y = C : (V • e*) . (5) 

If the transformation strains are known, then Eq. (4) can be solved for £ using 
Fourier transforms or Papkovich-Neuber potentials [31,32], 

Consider the particular case of an elliptic region in an infinite body, undergoing 
change of form that, without the constraint imposed by the surrounding ma- 
terial would result in an arbitrary homogeneous strain e t . The elliptic region, 
defined by 



has volume V and is bounded by the surface S. The fictitious body forces are 
in this case nil everywhere except on the surface S where they are 



n being the normal to the ellipsoid and 6 (•) the Dirac delta function. 

Within the ellipsoid, the total strain field £ = £ el + e t is uniform and can be 
expressed as [29] 

£ = P -.C-.y ( 8 ) 

where the fourth order tensor P results from solving Eq. (4) with b* from Eq. 
(7), and depends on the elastic properties of the material and geometry of the 
ellipsoid only. The derivation of the tensor P can be found in Refs. [31,32], 
The stress field is then obtained as 

a = C : £ el = C :(£-£*) = - (C - C : P : C) : £* = -Q : e* (9) 
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with 


Q = C - C : P : C. 


( 10 ) 


Eshelby’s inclusion problem Eshelby [30] showed that the eigenstrain 
problem can be used for the determination of the stress and strain fields due 
to an elliptic inclusion. Consider an elliptical inclusion with constitutive tensor 
C* in a homogeneous infinite solid with constitutive tensor C. Suppose next 
that the solid is loaded by a uniform stress or strain at infinity, a°° and £°°, 
respectively. The stress and strain in the solid can be expressed as 

<7 = <7°° + a and £ = £°° + £ (11) 

where a°° and £°° are the uniform stress and strain tensors induced in the 
solid if the inclusion was not present, while a and £ represent a perturbation 
due to the presence of the inclusion. 

The perturbation due to the inclusion can be computed using the eigenstrain 
analogy. In fact, the stress held a = cr°° solves the equilibrium equations 
everywhere in the solid, except in the inclusion, where the error in the stress 
is uniform: (C* — C) : £°°. This suggests that the stress state can be corrected 
by using a transformation strain inside the inclusion. For the inclusion and 
the transformed region to be equivalent, the stress in both cases must be the 
same. For the case of the hypothetical transformed region, the stress is 

a = C : £ el = C : (e - e*) (12) 

while in the case of the inclusion, the stress tensor is obtained as 

a = C* :£. (13) 

For the transformed region to be equivalent to the inclusion, then Eqs. (12) 
and (13) can be equated, resulting in 

C : (e - e*) = C* : e. (14) 

Decomposing the strain in its two components £°° and £, and considering Eq. 
(8), then Eq. (14) is obtained as 

(C* - C) : £°° = [C - (C - C*) : P : C] : e*. (15) 

For a void, C* = 0 and Eq. (15) reduces to 

e % = Q 1 : cr°°. (16) 
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The determination of £* in Eq. (16) is an important result, since the strain 
and stress at the cavity wall follow as 


£ = £°° + £ = £°° + P : C : £* = £°° + P : C : Q~ V°° (17) 

a — C : £. (18) 


The interaction energy can be defined as the energy released by the introduc- 
tion of a cavity [29]. Following from Eshelby’s work for a linear material [30], 
the interaction energy for a solid loaded at infinity can be written as 


Eini 


lv 



dV 


re* 

V I <7°° : de l 
Jo 


(19) 


where V is the volume of the cavity. 


The solution for a crack was obtained from the solution for an ellipsoidal 
cavity by Laws [29] . Laws considered first an infinite elliptic cylinder by letting 
c — > oo, then expressed the interaction energy per unit length of the cylinder 
in a form similar to 


Ej„t = tt are 


re 1 

/ cr°° : efe*, with e = 
Jo 


(20) 


and proceeded to a crack by making e — > 0. Since the tensor Q becomes 
singular when e — » 0, but not the product eQ _1 , some care has to be taken. 
Eq. (20) can be transformed in 


Em = vr a 2 / <7°° : eQ 1 : da c 


>0 


( 21 ) 


Making e — » 0 and defining A = eQ , the interaction energy can be expressed 
as 


EJi 77 + — 7 r a 


(T°° : A :da c 


10 


(22) 


where the non-zero components of the tensor A were calculated by Laws [29] . 
Assuming a non-linear shear behaviour, Eq. (22) can then be written as 


Eint — ( AqoCTqo + A o S Too + 2 


. o 2 


22 u 22 


v 23 ' 23 


T 12 A. 0 12 dT 12 


— 2 7ra ' ^ 


r 712 


^ 12^712 


2 na 


E° 22 o 22 + A ° 3 t 23 + x (712 ) 


(23) 


where x ( 712 ) is defined as 


X ( 712 ) = 2 



' 72^712 


(24) 



and, from Laws [29], 


A 


O 

22 



U 21 \ 

E\\ J ' 


(25) 


It is important to take into account some simplifications of the previously 
presented formulation: 

(i) Laws [29] derivation of the crack tensor A is only approximate. 

(ii) The superposition principle has been used, and it is not strictly valid for a 
non-linear shear behaviour. It is not clear if an expression for the interaction 
energy such as Eq. (23) could also be derived without using superposition. 


In-plane shear contribution to the critical energy release rates The 

in-plane shear contribution to the critical energy release rate can be com- 
puted from Eq. (23) with a- 22 = t 2 3 = 0. Suppose that the dimension of the 
crack in the longitudinal direction is o^. Proceeding like Dvorak and Laws 
[9], the change in cracked area for crack propagation in the transverse direc- 
tion is dA = 2az,<9a 0 and the critical energy release rate for the case of crack 
propagation in the transverse direction is obtained as 


rp da,LEi nt 

G ° = dA 


1 d Eint 

2 da Q 



(26) 


where 7 “ 2 is the engineering shear strain at failure. For crack propagation in 
the longitudinal direction, the change in cracked area is dA = 2a Q daL and the 
critical energy release rate is defined as 


da L E int 

Gc = dA 


E„ 


int 


2 CL 0 



(27) 


3.1.2 In- situ effect 

A failure criterion to predict matrix cracking under the presence of both in- 
plane shear and transverse tensile stresses should represent the ‘in-situ’ effect 
occurring in laminated composites. The in-situ effect, originally detected in 
Parvizi’s tensile tests of cross-ply glass fiber reinforced plastics [33], is char- 
acterized by higher transverse tensile and shear strengths of a ply when it 
is constrained by plies with different fiber orientations in a laminate, when 
compared with the strength of the same ply in a unidirectional laminate. The 
in-situ strength also depends on the number of plies clustered together, and 
on the fiber orientation of the constraining plies. 

The orientation of the constraining plies and the number of plies clustered 
together also affect the crack density and the stiffness reduction of the cracked 


9 




Fig. 1. Transverse tensile strength as a function of number of plies clustered together, 
with models from Dvorak [9] based on experimental data from Wang [34] . 

ply. Wang’s [34] tests of (0/90 n /0) (n = 1,2, 3, 4) carbon/epoxy laminates 
have shown higher crack densities for thinner 90° layers. The reduction of the 
elastic properties of a cracked ply is normally predicted using elastic analyses 
of cracked plies [15,35] or Continuum Damage Models [36-39]. 

The in-situ effect is illustrated in Fig. 1, where the relation between the in-situ 
transverse tensile strength and the total thickness of the 90° plies clustered 
together is represented. 

Accurate in-situ strengths are necessary for any stress-based failure crite- 
rion for matrix cracking in constrained plies. Both experimental [34,40,41] 
and analytical methods [9,35,42] have been proposed to determine the in-situ 
strengths. In the following, the in-situ strengths are calculated using fracture 
mechanics solutions for the propagation of cracks in a constrained ply. 


3.1.3 Fracture Mechanics Analysis of a Cracked Ply 

The failure criterion for predicting matrix cracking in a ply subjected to in- 
plane shear and transverse tension proposed here is based on the fracture 
mechanics analysis of a slit crack in a ply, as proposed by Dvorak and Laws 
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Fig. 2. Slit crack geometry (after Dvorak [9]) 

[9]. The slit crack represents a manufacturing defect that is idealized as lying 
on the (1,3) plane, as represented in Fig. 2. It has a length 2 a Q across the 
thickness of a ply, t. Physically, this crack represents a distribution of matrix- 
fiber debonds that may be present in a ply as a consequence of manufacturing 
defects or from residual thermal stresses resulting from the different coefficients 
of thermal expansion of the fibers and of the matrix. Therefore, the slit crack is 
an ‘effective crack,’ representing the macroscopic effect of matrix-fiber debonds 
that occur at the micromechanical level [34], 

The transverse tensile stress 022 is associated with mode I loading, whereas 
the in-plane and transverse shear stresses T12 and T23 respectively are associ- 
ated with mode II loading. The crack represented in Fig. 2 can grow in the 
1 (longitudinal, L) direction, in the 3 (transverse, T) direction, or in both 
directions. 

The components of the energy release rate for the crack geometry represented 
in Fig. 2 were determined by Dvorak and Laws [9] for a linear orthotropic 
material and an extension of their analysis for non-linear shear behaviour has 
been presented here. For mixed-mode loading, the energy release rate for crack 
growth in the T and L directions, G(T ) and G(L), respectively, are given by 

G(T) = ^ [vj^22 (T 22 + ViA 23 T 23 + TlIlX (712)] , . 

G(L) = K?a> 1 2 + &a; 3 t | 3 + e,nx (712)] 

where it can be observed that the energy release rate G(L) for longitudinal 
propagation is a function of the transverse slit size and that it is not a function 
of the slit length in the longitudinal direction. 

The parameters rji, i = I, II, III in Eq. (28) are stress intensity reduction 
coefficients for propagation in the transverse direction, and the parameters 
£i, % = I, II, I II are reduction coefficients for propagation in the longitudi- 
nal direction. These coefficients account for the constraining effects of the 
adjoining layers on crack propagation: the coefficients are nearly equal to 
1.0 when 2 a Q -C t, and are less than 1.0 when a Q ~ t. Experimental re- 
sults [41] have shown an increase in the in-situ transverse tensile strength 
of [i^/OOfijs, 9 = 0°,30°,60°, laminates for increasing stiffness of adjoining 
sublaminates ± 9 . This implies that the value of the parameter 7 ^ decreases 
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with increasing stiffness of adjoining sublaminates. Considering that a trans- 
verse crack can promote delamination between the plies, Dvorak and Laws 
[9] suggested that the effective value of rj l can be larger than obtained from 
the analysis of cracks terminating at the interface, and suggested the use of 
Vi = & = !• 

The mode II and III components of the energy release rate are combined in a 
shear mode, Gsh , as 

Gsh = Gn + Gm • (29) 

Such an approach was initially proposed by Li and Sen [43] and Li [44], and 
used in the simulation of delamination using the Virtual crack Closure Tech- 
nique (VCCT) and decohesion finite elements [2], By combining modes II and 
III, it is not necessary to track the relative orientation of the crack front with 
respect of the in-plane displacement jumps. In addition, no conclusive evi- 
dence is available showing that Gu c and Gjjj c are different. In fact, there is 
no standard for measuring Gm c . Furthermore, there is no mixed-mode test 
method for mixed mode II and mode III loading. 

The components of the energy release rate are then obtained for the In- 
direction using Eq. (28) with r/,; = 1: 

Gi(T) = ^A> 2 2 (30) 

jr n r 1 

Gsh(T) = ~Y [A 23 t 2 2 3 + x (712)] . (31) 

The corresponding components of the fracture toughness are given as 

Gic(T) = ^ (Fj) 2 (32) 

GWO = (33) 

where Y? is the in-situ transverse tensile strength, and 7 ]' 2 is is the in-situ 
in-plane shear ultimate strain. 

For propagation in the longitudinal direction, the mode I and mode II com- 
ponents of the energy release rate are 

Gi(L) = — ^A 22 ct 22 (34) 

•77-/7 r 1 

Gsh(L) = Y [ A 2s4 + X ( 712 )] (35) 

and the components of the fracture toughness are 

G, C (L) = ^A° 22 (y;) 2 (36) 
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( 37 ) 


GsHc(L) = — pX (71*2, is) ■ 


Having obtained expressions for the components of the energy release rate and 
fracture toughness, a failure criterion can be applied to predict the propaga- 
tion of the slit crack represented in Fig. 2. Under the presence of in-plane and 
transverse shear, as well as transverse tension, the critical energy release rate 
G c depends on the combined effect of all microscopic energy absorbing mech- 
anisms such as the creation of new fracture surface. Relying on microscopic 
examinations of the fracture surface, Hahn [45] observed that the fracture 
surface topography strongly depends on the type of loading. With increasing 
proportion of the stress intensity factor A'//, more hackles are observed in 
the matrix, thereby indicating more energy absorption associated with crack 
extension. Hahn proposed a mixed-mode criterion written as a first-order poly- 
nomial in the stress intensity factors Kj and Kjj. Written in terms of the mode 
1 and mode II energy release rates, the Hahn criterion is 


(1-9) 


\ 


Gj(i) , G r (i) , G S h(i) 


G, 


+ 9 


n. 


+ 


= 1, i = T,L 


(38) 


where the material constant g is defined from Eqs. (36) and (37) as g — 
Gic/Gshc which leads to the following expression for g: 


a ; 2 (Q 2 


(39) 


A failure index for matrix tension can be expressed in terms of the ply stresses 
and in-situ strengths by substituting either Eqs. (30)-(33) or (34)-(37) into the 
criterion in Eq. (38) to get 


LaRC04 #1 


FI m = 


(1-9) 


&22 , 

YT+9 
1 is 



A 23 r l 3 + X (712) 

X (7l2,w) 


1. (40) 


The criterion presented in Eq. (40), with linear and quadratic terms in <722, a 
quadratic term in T23 and a term on the in-plane shear internal energy, y (712), 
is similar to the criteria proposed by Hahn [45] , Liu [24] (for transverse tension 
and in-plane shear), and Puck [6]. It can be observed that using g — 1 in Eq. 
(38) results in the linear version of the criterion proposed by Wu and Reuter 
[46] for the propagation of delamination in laminated composites 


Gj G n 
G Ic G IIc 


1 . 


(41) 


Furthermore, using g — 1, assuming linear in-plane shear, and neglecting 
r 2 3, Eq. (40) reverts to the well-known Hashin criterion [11] for transverse 


13 




Fig. 3. Unidirectional specimen with initial crack (after Dvorak [9]) 

matrix cracking under both in-plane shear and transverse tension, where the 
ply strengths are replaced by the in-situ strengths 


FI m = 



(42) 


Finally, the non-linear term in Eq. (40) is also found to be similar to the 
strain-energy based criterion proposed by Sandlin [47], later used by Chang 
and Scott [48]. 


3.1.4 Application to unidirectional laminates 

The application of the fracture mechanics analysis of a cracked ply to unidirec- 
tional laminates yields expressions relating toughness values to crack dimen- 
sions which can subsequently be used to relate the strength of thick embedded 
plies to that of a unidirectional laminate. 

Dvorak and Laws [9] regarded the fracture of a unidirectional specimen as 
the fracture of an unconstrained thick ply, with the critical initial slit crack 
located at the surface of the laminate. For tensile loading, the crack can be 
located at the edge of the laminate, which increases the energy release rate 
when compared with a central crack. In the case of shear loading, there is no 
free edge, so the crack is a central crack, as shown in Fig. 3. The defect size 
is 2 a Q and is considered to be much smaller than the ply thickness, 2 a Q -C t. 

For unidirectional laminates, the crack will grow unstably in the transverse 
direction [9], and Eqs. (32), (33), (36) and (37) apply with a geometric fac- 
tor which is obtained from the classic solution of the free edge crack [5,49], 
resulting in 

G Ic (T) = 1.12 2 7ro 0 A°(F T ) 2 (43) 

Gshc{T)=™ oX {'Y? 2 ) (44) 

where Y T is the transverse tensile strength and Y 12 i s the in-plane shear strain 
at failure, for a unidirectional laminate. 

The toughness ratio g can also be calculated in terms of the unidirectional 
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Fig. 4. Geometry of slit crack in a thick embedded ply subjected to tension and 
shear loads 

properties by using Eq. (43) and (44), resulting in 


9 = 


Gic 

GsHc 



(45) 


3.1.5 In- situ strength of thick embedded plies 

A thick ply is defined as one in which the length of the slit crack is much 
smaller than the ply thickness, 2 a Q <C t, as illustrated in Fig. 4. The minimum 
thickness for a thick ply depends on the material used. For E-glass/epoxy 
and carbon/epoxy laminates, Dvorak and Laws [9] calculated the transition 
thickness between a thin and a thick ply to be approximately 0.7mm , or about 
5 to 6 plies. 


For the geometry represented in Fig. 4, the crack can grow in the transverse or 
in the longitudinal direction. Comparing Eqs. (30) and (31) to Eqs. (34) and 
(35), however, indicates that the energy release rate for the crack slit is twice 
as large in the transverse direction as it is in the longitudinal direction. Since 
Eqs. (30) and (31) also indicate that the energy release rate is proportional to 
the crack length, the crack will grow unstably in the transverse direction. Once 
the crack reaches the constraining plies, it can propagate in the longitudinal 
direction, as well as induce a delamination. 


Crack propagation is predicted using Eq. (40), and the in-situ strengths can 
be calculated from the corresponding fracture toughness, as expressed in Eqs. 
(32) and (33). For mode I, the in-situ transverse tensile strength is defined 
from Eq. (32) as 


Y: = 


2 GUT) 

\ na 0 A 


O 

22 


(46) 


and, taking into account Eq. (43), 


Y't = 1.12V2Y T . 


(47) 


For a thick embedded ply loaded in pure in-plane shear, the expressions get 
more complex due to shear non-linearity. The in-situ in-plane ultimate strain 
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is obtained as 

7i2,* s = X I 

and, taking into account Eq. (44), 

7l2i is = X 


2GWT)\ 

(48) 

TICLq J 

[2x(7i" 2 )]. 

(49) 


Considering the constitutive law for the shear behavior to be expressed by the 
function fc l such that 

r = fcL ( 7 ) , (50) 

then the in-situ in-plane shear strength is defined as 

St = f cl (tEm,) ■ (51) 

As a particular case of non-linear in-plane shear behaviour, consider the poly- 
nomial relation between the shear strain and the shear stress proposed by 
Hahn and Tsai [50], 

712 = 7^072 + /foi 2 - ( 52 ) 

Ct 12 

With the non-linear shear law expressed in Eq. (52), Eq. (48) reverts to 


GsHciT) 




and Eq. (44) can be written as 


GsHciT) = 2na 0 



(53) 


(54) 


Equating Eqs. (53) and (54), the in-situ shear strength of a thick embedded 
ply, S^ s , can be related to the shear strength of a unidirectional laminate: 



(55) 


The in-situ shear strength of an thick embedded ply, S^ s , is the positive, real 
root of Eq. (55): 


. V 

'1+/3I 

(X ,2 + 18/3(S0 4 ) 

1 ( G 12 ) 2 - 1 

\ 3 PG 12 


(56) 


It can be observed from Eqs. (47) and (56) that the in-situ strengths of thick 
embedded plies — Y?, and S -' s — are independent of the ply thickness, as has 
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Fig. 5. Geometry of slit crack in a thin embedded ply 

been observed by Dvorak and Laws [9] and Leguillon [51], and as was shown 
in Fig. 1. The general expression for Sh in Eq. (56) can be written for a linear 
material by letting f3 tend to zero, in which case the in-plane shear strength 
come as 

S^ s = v / 2S' i for a linear shear law. (57) 

Eq. (57) is equal to the one obtained by Dvorak and Laws [9] and Davila et 
al. [5] for a linear shear behaviour. 


3.1.6 In-situ strengths of thin embedded plies 


Thin plies are defined as having a thickness smaller than the typical defect, t < 
2 a 0 , so the slit crack represented in Fig. 2 extends across the entire thickness 
t of the ply, as represented in Fig. 5. 


In the case of thin plies, crack defects can only grow in the longitudinal (L) 
direction, or trigger a delamination between the plies. The in-situ strengths can 
be calculated from the components of the fracture toughness, as expressed in 
Eqs. (36) and (37). For transverse tensile loading, the corresponding strength 
is expressed as 

yT _ 8G Jc (L) 

“ \ JrtAfc ' 

For a thin embedded ply loaded in pure in-plane shear, the ultimate in-plane 
shear strain is defined as 



7l2iis — X 


8 Gshc(L) 


Tit 


and the in-situ in-plane shear strength is defined as 

St = fcL (7?2,i.) • 


(59) 


(60) 


Considering again Eq. (52) as a particular case of non-linear in-plane shear 
behaviour, Eq. (59) can be written as 



Gshc(L) 

Tit 


(61) 


The in-situ shear strength of a thin ply, Sf s , is the positive real root of Eq. 
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( 61 ): 


St = 


l i+p A8G SHc (L) {G ) 2_ 1 

7 rt 


(62) 


\ 3f3G\2 

It can be observed from Eqs. (58) and (62) that the in-situ strengths are 
dependent on the thickness t. 


Eq. (62) can be written for a linear material by letting (3 tend to zero, in which 
case the shear strength come as 

Sf s = \J ^ ^ 1< iy ^ for a l' near shear law. (63) 

Eq. (63) is equal to the one obtained by Dvorak and Laws [9] and Davila et 
al. [5] for a linear shear behaviour. 

If in the absence of specific data, the toughness values Gi c (L) and Gshc(L) 
can be assumed to have the values measured by standard Fracture Mechanics 
tests, such as the DCB for mode I and the ENF test for mode II. Using Eq. 
(58), Dvorak and Laws [9] obtained a good correlation between the predicted 
and experimentally obtained in-situ tensile strengths of both thick and thin 
90° plies in [0/90 n /0] laminates, as was shown in Fig. 1. 


3.2 Compressive matrix failure 


Matrix compression specimens fail by shear, which would suggest that the an- 
gle of the fracture surface with the through-the-thickness direction, as shown 
in Fig. 6(a), should be a 0 = 45°, i.e., along the plane of the maximum shear 
stress. However, experiments indicate that the angle of fracture under uniaxial 
compression is generally a Q = 53 ± 2° for most technical composite materi- 
als [6,52,53]. The fact that a 0 > 45° can be explained by the existence of 
a compressive stress acting on the potential fracture surfaces, and its associ- 
ated friction. The magnitude of the compressive stress, and hence the friction 
stress, is maximum for a fracture surface with a Q = 0° and decreases mono- 
tonically until a Q = 90°, in which case the compressive (and friction) stress is 
zero. Although the shear stress is maximum for a 0 = 45°, the friction stress 
which opposes fracture, decreases with larger values of the angle a 0 . As a 
result, fracture is expected for values of a 0 larger than 45°, where a critical 
combination of shear and normal stress acts. 

Experimental evidence on the fracture surface of specimens failing by matrix 
compression suggest that the Mohr-Coulomb (M-C) criterion is applicable to 
the matrix compression failure; Puck [13] and Puck and Schneider [14] were 
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the first to propose a matrix failure model for composites based on the M-C 
criterion. 


The M-C criterion is commonly used in applications where fracture under ten- 
sion loading is different from fracture under compression loading, such as in 
soil mechanics or in the fracture of cast iron. The application of the M-C crite- 
rion to multiaxial failure of epoxy resins was studied by Kawabata [54] based 
on correlation with his own test results. While studying the failure of chopped 
glass- fiber /epoxy mat laminates under confining pressures, Bochler and Raclin 
[55] found the Tsai-Wu criterion to be inadequate, and formulated a shearing 
criterion based on the M-C criterion that fit his experimental measurements 
well. Taliercio and Sagramoso [56] used the M-C criterion within a non-linear 
micromechanical model to predict the macroscopic strength properties of fiber 
composites. 


The M-C criterion is represented geometrically by the diagram illustrated in 
Fig. 6(b). The Mohr’s circle represents a state of uniaxial compression. The 
angle of the plane of fracture is a 0 . The M-C criterion postulates that in a 
state of biaxial normal stress, fracture occurs for any Mohr’s circle that is 
tangent to the M-C fracture line. 


DiLandro and Pegoraro [57] explain the role of internal friction on the strength 
of carbon-fiber composites by noting the absence of chemical bonds between 
fiber and matrix, and that adhesion is attributed to Van der Waal’s interac- 
tions. When subjected to an external load, the shear slipping of the two phases 
is prevented until the shear stress at the fiber-matrix interface reaches a lim- 
iting value. DiLandro and Pegoraro also note that the friction coefficient is an 
empirical factor that encompasses all chemical-physical interactions, includ- 
ing the thermal residual shrinkage of the matrix around the fiber. Larson and 
Miles [58] examined the relative effects of interfacial friction and roughness on 
the length of interfacial sliding which proceeds from the tip of an impinging 
fracture oriented perpendicular to the interface. According to Larson, sliding 
is key to the cracking behavior of fibrous brittle matrix composites in that it 
affects the stress concentration on the fibers, the matrix crack spacing, and, 
therefore, the global toughness of a composite material. 


For a general loading situation, shown in Fig. 6(c), the angle of the fracture 
plane with the through-the-thickness direction, denoted as a, might assume a 
different value than the one for pure compression (a 0 ). The orientation of the 
fracture plane depends on the particular combination of shear (r T and r L ) and 
normal (a n ) tractions for each value of a, shown in Fig. 6(d). In plane-stress 
formulations, the tractions at potential fracture planes are obtained from the 
in-plane shear stress T 12 , the compressive stress cr 2 2 , and the fracture plane 
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Fig. 6. (a) Pure transverse compression failure, for a CFRP specimen; (b) stresses 
in the fracture plane; (c) geometrical representation of the Mohr-Coulomb criterion; 
(d) fracture plane for a 3D stress state; 

angle a by using the transformation equations 

o n = °y (1 + cos ( 2 q 0) 

< r 1 ~”Dpsin(2 a) (64) 

t l = T12 cos (a) . 

In a 3D formulation, the tractions are obtained from the components of the 
stress tensor and the fracture plane angle a: 

O n = t ^ C ° S ( 2a ) + T 23 S111 ( 2q 0 

< r T = ~~ 2 — ~ s * n ^ 01 ^ T23 cos (65) 

r L = T12 cos (a) + r 31 sin (a) 


where a is comprised in the interval ] — 7T, 7t] . 


The M-C failure criterion is expressed in terms of the tractions in the frac- 
ture plane, and can be written in several forms. Considering first the case 
where t l = 0 (i.e. , T\2 = 0 for the 2D formulation), the M-C criterion can be 
expressed as 


r 


+ V T On 


= s 1 


( 66 ) 


where 7] T is a friction coefficient and S T is the transverse (to the fibers) shear 
strength. The geometrical representation of this criterion in a (cr,|r|) space is 
a line with negative slope (— r] T ), shown in Fig. 6(b). In this figure, the Mohr 
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circle corresponding to the case of pure compression is represented through 
a semi-circumference tangent to the M-C criterion’s line. The slope of the 
M-C criterion’s line can be related to the angle of the fracture plane in pure 
transverse compression, a 0 , through 

1 

tan(2o; 0 ) = (67) 

V 

Furthermore, writing Eq. (66) for a pure compression case establishes the 
relation between S T , Y c and a Q : 

S T = Y° cos (a 0 ) (sin (a G ) + cos ] . (68) 

\ tan (2a 0 ) J 

The angle a Q can be easily determined from simple compression tests. The 
parameters rf and S T are calculated from a Q using Eqs. (67) and (68). 


The M-C criterion (Eq. (66)) can be expressed in several forms, namely con- 
sidering that friction affects (increases) the strength, 


T 


= 1 , 


S T - if a, 

or that is affects (decreases) the applied stress, 

( r T + rfa n ) 


(69) 


S T 


= 1, 


(70) 


where the operator (•) is the Me Cauley bracket defined by 


0 4= x < 0 

(x) — < , igR. 

x 4= x > 0 


( 71 ) 


Motivated by those two different forms of expressing the same criterion for 
t l = 0, Puck and Schiirmann [6,52] initially proposed for the general case 

(r L f 0 ) 


FI m = 


T 


S T 


ifa n 


+ 


T 


S L 


r] L a n 


< 1 


(72) 


whereas it was proposed, first for the LaRC02 [8], and after for the LaRC03 1 
[5] failure criteria 


FI m = 



( r L + 'r] L a n ) ' 


S L 


< 1 


(73) 


1 For the LaRC03 criteria [5], the in-situ longitudinal shear strength SF was con- 
sidered in equation 73 instead of S L 
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where i] T (r n and r] L a n are shear stresses due to friction. While Puck and 
Schurmann (Eq. (72)) consider that the compression stress (a n ) increases the 
effective strength, LaRC02 and LaRC03 (Eq. (73)) consider that the compres- 
sion stress reduces the effective shear stress. For reasons of simplicity, Puck 
and Schurmann [6,52] transformed Eq. (72) into 


FI m = 




+ 


(W 


(S T ) 2 - 2 ifS T a n ( S L f - 2rj L S L a, 


< 1. 


(74) 


Eq. (74) is obtained by expanding the denominators in Eq. (72) and neglecting 
the terms in (cr n ) 2 . Puck and Schurmann [6] have shown that it is possible to 
solve Eq. (74) in a closed form for the angle of the fracture plane a in the case 
of a plane stress situation. 

To obtain the friction coefficient rj L in the absence of experimental data, Puck 
and Schurmann [6] suggest the relation 


r r p 

Tj T] 

~S l = ~S f 


(75) 


where r] T is computed using Eq. 67 and S T is computed using Eq. 68. This 
suggestion has also been considered for the LaRC02 [8] and LaRC03 criteria 
[5], 

The comparison of the criteria expressed in Eqs. (72), (73) and (74) is shown in 
Fig. 7 for a E-glass/DY063 epoxy used in the WWFE, with material properties 
presented in Table 2. Also shown in Fig. 7 is the much simpler Sun [23] criterion 


FI m = 



+ 


Tl2 


S L 


V L V22, 


< 1. 


(76) 


The Sun criterion also represents a similar trend in the failure envelope, al- 
though it is not based on a physical mechanism of failure. 

While the initial and final Puck criteria (Eqs. (72) and (74)) yield similar 
results, LaRC03 criterion (Eq. (73)) is less conservative. This is related to the 
fact that the effect of friction is over-estimated in Eq. (73). Indeed, affecting 
the shear tractions by a friction term as in Eq. (73) over-estimates the friction 
forces whenever both r T and r L are acting simultaneously [59]. As Fig. 8(a) 
represents, supposing a very simple case with isotropic friction (■ r] T = r] L ), the 
friction stresses are over-estimated by a factor of \/2 when using Eq. (73). 

It is interesting to notice the effect that an orthotropic friction model has on 
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Table 2 

Mechanical properties of a E-glass/DY063 epoxy used in the WWFE, from Soden 
et al. [60] 


G \ 2 

S L 

X c 

Y c 

Oi 0 

(GPa) 

(MPa) 

(MPa) 

(MPa) 

(°) 

5.83 

73 

800 

145 

53 



Fig. 7. Failure envelopes for transverse compression and inplane shear 


* T fnc =/7 r x<T„xcos(0) 

*fric = ? /XCX„ 

with T] = yj \rf cos(#))~ + (r] L sin (6*))“ 

^Tfnc = 7 L x<r B xsin(0) 

(b) 

Fig. 8. (a) Overestimation of the friction stress; (b) model that does not overestimate 
the friction stresses 



LaRC03 criterion (Eq. (73)). A reasonable model for orthotropic friction is 


fric. 

r T 

fric. 

n 




r] T 0 
0 r] L 



(77) 


23 


x 12 (MPa) 



120 

100 

80 

60 

40 

20 

0 


-150 -130 -110 -90 -70 -50 -30 -10 

0 22 (MPa) 


Fig. 9. Failure envelope for the LaRC04 matrix failure criterion, using Eq. 81 and 
several values for the angle of the fracture plane, a 


or, simplifying, 


7j/"' = a n r] T cos ( 9 ) 
-/■ = tr „n L sill (0) 


(78) 


where 9 is the angle formed by the shear-traction vector r and the transverse 
direction in the fracture plane, i.e. , 9 = arctan (y L /r T ^j. 

The LaRC03 criterion given by Eq. (73) can be modified using Eqs. (78) to 
account more accurately for the friction terms, which results in 


FI m = 


(L 

t t 

+ r) T O n COS {9)y 

2 

+ 

« 

r L 

+ r/ L o n sin (#)) ^ 

\ S T J 

i S L / 


< 1 


(79) 


Interestingly, the modified criterion given by Eqs. (79) yields an envelope that 
is coincident with Puck’s initial criterion (Eq. (72)), as shown in Fig. 7. 

The use of one of Eqs. (72), (73), (74), and (79) for the failure criterion implies 
the use of the set of Eqs. (65) for the transformation of stresses. In turn, 
this means that the fracture angle has to be known. In theory, it could be 
analytically determined as 


a : 


OFIm 

da 


(80) 


Puck and Schiirmann [6] have solved Eq. 80 for the criterion in Eq. (74) in 
a plane-stress situation, but it is not possible to do the same with Eqs. (72), 
(73) and (79) for 3D stress states. However, Fig. 9, drawn using the data from 
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Table 2 and Puck’s initial criterion (Eq. (72)), shows that it is possible to 
draw the envelope with reasonable accuracy by using a very small number of 
trial angles. 

Since it is not possible to obtain a closed-form solution for the fracture angle in 
a generic 3D situation, there does not seem to be a reason to drop the physical 
significance of the criteria expressed through Eqs. (72), (73) and (79). As it 
was demonstrated, Eq. (73) over-estimates the friction stresses. Since Eqs. 
(72), (79) yield similar results, but Eq. (72) is simpler, the later is proposed 
for LaRC04, 


LaRC04#2 FI U = f ^ ) + ( A — L ) <1 (81) 

\S T — r) T a n J \S{ s — rj L a n J 

where it should be noted that the in-situ longitudinal shear strength S^ s is 
used, as was proposed for LaRC03 [5]. 


4 LaRC04 criteria for fiber failure 


4- 1 Tensile fiber failure 


The LaRC04 criterion for fiber tensile failure is a non-interacting maximum 
allowable stress criterion. Consequently, the LaRC04 failure index for fiber 
tensile failure is simply 

LaRC04 #3 FI f = ^ < 1. (82) 


4-2 Compressive fiber failure 


Fiber compression failure is a field where significant research is still being per- 
formed. For matrix compression failure, a relatively simple mechanical model 
as the one proposed by Puck and Schurmann [6,13,14,52] seems to correctly 
represent failure, and can be easily expressed as a failure criterion that can 
be incorporated in numerical codes. Depending on the material, different fiber 
compressive failure modes are possible [61]: 

Microbuckling. This failure mode consists of the microbuckling of the fibers 
in the elastic matrix. The first mechanical model for this failure mode can be 
tracked back to Rosen’s work [62] where the fibers are represented by infinite 
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beams 2 in an elastic matrix and failure is attained when the compressive load 
equals the buckling load. This model provides an upper bound for the failure 
stress, as it generally predicts a failure stress typically two to three times larger 
than the experimental one (for carbon reinforced composites [61]). Models 
based on microbuckling have been widely studied over the last decades. For 
these models, the matrix shear properties as well as material imperfections 
play an important role. 

Kinking. Kinking can be defined as the localized shear deformation of the 
matrix, along a band. Typically, the fibers break at the edges of the band, and 
sometimes also in the interior. It should be noted that some authors consider 
kinking as a consequence of microbuckling, while others consider it as a sepa- 
rate failure mode [61]. Argon [7] was the first to develop a mechanical model 
for fiber kinking as a separate failure mode. For Argon, failure is the result 
of matrix shear failure, prompted by an initial fiber misalignment. For this 
model, and those that follow it, matrix elastic behaviour and initial material 
imperfections play an important role. 

Fiber failure. Fiber failure can be expected to occur for fibers with low 
compressive strength, such as Ararnid, but is not expected to happen for 
carbon, glass or boron fibers [61]. 

Rosen [62] was the first to propose a mechanical model to describe fiber mi- 
crobuckling. By minimizing the internal energy, Rosen obtained the fiber’s 
buckling stresses. The analysis was performed for two instability modes: (i) 
an extension mode, where the deformed shapes of adjacent libers are in oppo- 
sition of phase; and (ii) a shear mode, where the deformed shapes of adjacent 
fibers are in phase. Rosen found that for composites with a high fiber volume 
fraction, the shear mode is critical, and the associated stress is 

X c = T%- (83) 

1 ~ V f 

where G m is the shear modulus of the matrix and Vf is the fiber volume frac- 
tion. Considering the rule of mixtures, Eq. (83) reduces to X c = G 12 , he., the 
shear modulus of the composite. This relation (Eq. (83)) was later modified to 
account for an elastic-perfectly plastic resin [63]. In fact, several modifications 
were attempted, in order to incorporate less restrictive hypotheses (see Ref. 
[61] for a state of the art review). However, Rosen’s approach yields smaller 
failure stresses than similar models assuming linear elasticity and straight 
fibers [61] and is simpler. Still, there is little success in predicting the fail- 
ure stress of advanced composites using Rosen’s result. The problem is that, 
when compared to experimental data, Eq. (83) gives results typically 1.5 times 

2 In fact, Rosen’s approach is 2D and the fibers are thus represented by layers 
(plates), and not beams 
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higher for Boron composites, 2 to 3 times higher for carbon composites and 4 
times higher for glass composites [61]. It was with the introduction of geomet- 
ric non-linearity and initial fiber misalignment that the prediction got closer 
to the experimental results. 

Schultheisz and Waas [61] pointed out that most buckling models tend to repli- 
cate the model-composite that was studied by Greszczuk [64-67]. Greszczuk 
performed a series of experiments on model-composites whose reinforcements 
consisted of either (i) rods (of steel or aluminum) having diameters in the 
range 0.5 to 3.2mm ; or (ii) aluminum plates with thickness in the range 
of 0.3 to 1mm . The basic idea was to duplicate the 2D geometries used in 
the analytical models. Greszczuk found that the compressive strength of his 
model-composites exceeded Rosen’s prediction. However, when including the 
energy associated with bending of the fibers, Greszczuk found good agree- 
ment. He also concluded that while his model-composites with low-modulus 
matrix failed by microbuckling, those with intermediate-modulus matrix failed 
by longitudinal cracking (matrix cracking), and composites with high- modulus 
matrix failed through compression of the fibers. 

Most buckling models follow the 2D approach of Rosen [62], However, real 
technical composites are 3D structures. It has been suggested [68,69] that 
3D effects may be a cause of the reduction in the failure stress from Rosen 
model. Indeed, evidence of the 3D aspect of fiber microbuckling has been re- 
ported in the literature [70,71]. One important 3D consideration, discussed 
by Schultheisz and Waas [61], is the arrangement of fibers and matrix, which 
induce different types of interaction (among the fibers and matrix) during the 
buckling for different packing densities. Furthermore, the laminated construc- 
tion may lead to different fiber arrangements within the plane of a lamina 
and in the through-the-thickness direction. Fiber misalignment angles have 
first been reported to be smaller in the through-the-thickness direction [72], 
but more recent results suggest they are similar in magnitude [73]. Models in- 
corporating 3D aspects were attempted [67,74-78], and Schultheisz and Waas 
[61] concluded that FE analysis would be an excellent candidate to handle the 
3D aspect of microbuckling. 

In most high fiber-volume-fraction advanced composite materials, compressive 
failure is seen as a failure mode which is localized in a band across the speci- 
men in which the fibers have rotated by a large amount, and the matrix has 
undergone large shearing deformation — kinking. A schematic representation of 
a kink band is shown in Fig. 10(a). In the literature, kinking is often seen as a 
consequence of microbuckling, and not a failure mode itself. However, kinking 
is other times identified as an independent failure mode. On the discussion on 
whether or not kinking is a consequence of microbuckling, the main argument 
has to do with the orientation of the kink-band boundary. Indeed, if kinking is 
a consequence of microbuckling, then one would expect the kink-band bound- 
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(a) (b) 


Fig. 10. (a) Kinking band; (b) fiber misalignment frame 

ary to lie normal to the loading axis (original fiber direction), i.e., to lie in the 
plane of highest bending stresses. In Fig. 10(a), this would mean an angle f3 
equal to zero. However, it is found that in most cases / 3 lie in the range of 30°. 
On the other hand, the similarity between the kink bands and shear bands 
may suggest that shear is the main factor at the onset of kink-band formation. 
In this case, it would be expected that kink bands would occur in the planes 
of maximum shear stress, i.e., for (3 = 45°. Some thick kink bands have been 
found near 45° [61], but this is generally not the case. 

Another argument supporting kinking as a separate failure mode was intro- 
duced by Chaplin [79], who noted that microbuckling should occur everywhere 
in the composite at about the same time, whereas kinking in his experiments 
was clearly initiated from some kind of defect. 

Effendi [80] carried a set of experimental tests on different carbon-fiber com- 
posites and also carried an analytical buckling analysis, which included initial 
fiber waviness, and computed the stresses in the fiber and matrix constituents. 
He found that before the buckling load was attained, the stress levels in the 
constituents exceeded the respective failure stress. To refine the modelling, he 
conducted a numerical FE analysis, where a non-linear matrix behaviour was 
incorporated. The results confirmed that constituent failure happened first. 
Composites with small initial imperfections or weaker fiber would fail by fiber 
failure, while in composites with large initial imperfections or stronger fibers, 
the matrix would fail first. 

On their review of the state of the art on compression failure of composites, 
Schultheisz and Waas [61] concluded that the experimental observations on 
the kinking phenomenon support the contention that kink bands seen in high 


fiber-volume-fraction advanced composite materials occur via a mechanism 
that is different from the global microbuckling mode suggested by Rosen [62] . 
Furthermore, in advanced composites, kinking seems to be initiated by local 
microstructural defects, such as fiber misalignments and longitudinal (matrix 
or interfacial) cracking. Therefore, kinking is better understood by treating 
the problem of stress redistribution, including dynamics, and including both 
non-linear geometry and material response. 


Whether kinking is a result of fiber microbuckling or a separate failure mode, it 
is the most common failure feature observed after testing. The kink-band angle 
and kink width were studied by Hahn [81] for carbon fiber composites (CFRP), 
glass fiber composites (GFRP) and Aramid fiber composites (AFRC). The 
kink-band angle /3, and the band width w , were found to be the smaller for 
CFRP at room temperature (3 ~ 20° and w ~ 0.07 to 0.2mm). For GFRP 
no clear kink bands were observed at room temperature. However, at 100°C , 
a kink angle f3 ~ 30° and a width w ~ 1.2mm were observed. For AFRP 
tested at room temperature, the kink angle was found as (3 ~ 40° and width 
w ~ 0.45mm . Chaplin [79] noted the angle of rotation of the fibers in the kink 
band was twice the angle /?, so that no volumetric changes happened in the 
kinked region. 


Argon [7] was the first researcher to provide a model for kink-band formation. 
In his model, it is assumed that an initial fiber misalignment exists, which 
leads to shearing stresses between the fibers. The shearing stresses would act 
as to further rotate the fibers, which would in turn lead to further increase 
in the shear stresses. This closed-loop effect would then lead to failure. The 
main result from his analysis is the relation between the compressive failure 
stress, X c , the matrix in-plane shear failure stress, S L , and the initial fiber 
misalignment angle (f°: 


X c 


SP_ 

( p° 


(84) 


From Argon’s analysis, a kink-band angle (3 = 45° should also result. Several 
authors since reported the sensitivity of the compressive failure stress to the 
shear failure stress [82-85] . Budiansky [84] later extended Argon’s analysis to 




S L 

( p° + 


(85) 


where 7 “ is the shear strain at failure. Failure occurs when the shear failure 
stress is reached. 


4-2.1 2D Kinking model 

Consider a unidirectional composite with a misaligned region being com- 
pressed, as depicted in Fig. 10(b). The stresses in the misalignment frame 
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are 


<711+^22 <7 H — CT22 /q \ | • / 0 \ 

CTimim = 1 cos (2ip) + T 12 sin (2^) 

cr 2 m 2 m = cos (2 ip) - Ti 2 sin (2 (f) 

T lm2 m = — — sin (2 tp) + t \2 cos (2 f) . 


( 86 ) 


For failure under pure compression (an = —X° and a 2 2 = Ti 2 = 0), Eqs. (86) 
lead to 

i 

af ml m = —X c cos 2 (f ) , 

a 2 m 2m = X c sin 2 (f) and (87) 

T{m 2 m = X c sin (f) cos (f) . 


The stress state in the misalignment coordinate frame can now be used to 
check for fiber kinking. For a material with linear shear behaviour, replacing 
it in an appropriated matrix failure criterion leads directly to the expression for 
the specific value of the misalignment angle f at failure for a pure compression 
case — f c . For a material with non-linear shear behaviour, it will become clear 
at the end of this section that kinking can result either from (i) matrix failure 
(i.e. the verification of a matrix failure criterion), or (ii) instability, due to the 
loss of (shear) stiffness for larger shear strain values. 


Without loss of generality, suppose first the case of failure by the verification 
of an appropriate matrix failure criterion. 


Case 1 Kinking for pure compression as the result of matrix failure 


Using LaRC04 matrix compression failure criterion, Eq. (81), the expression 
for the specific value of the misalignment angle f at failure for a pure com- 
pression case (f c ) is obtained as 

(sin (f c ) cos (f c ) - r] L sin 2 (^ c )) = S L . (88) 

The angle f c is the sum of an initial misalignment with the rotation due to 
loading. Eq. (88) can be solved for f c [5], resulting in 


f c = arctan 


1 

T— 1 

1 

t-H 

-4 

(ft + V L ) 

^ |o 

l 

2 

( w + v L ) 

) 


(89) 


It can be observed that by neglecting rj L , and assuming f c to be small and 
constant, Eq. (88) yields Argon’s equation (Eq. (84)). In fact, assuming all 
the above but not that f c is constant yields Budiansky’s Eq. (85). 


From the constitutive law, the shear stress in the misaligned frame is obtained 
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as a function of the shear strain 


T lm2 m — fcL (7l m 2 m ) 


(90) 


and from the transformation equations (Eqs. (86)), the shear stress in the 
same frame for pure axial compressive failure is 

r 1 C m 2 m = ^ sin (2<^ c ) X c . (91) 

Equating Eqs. (90) and (91), the shear strain at failure for a pure axial com- 
pression case, 7im 2 m, is obtained as 

Yim 2 m = f cl Q sin (2 ip c ) X c ) . (92) 

For instance, for a material with a linear shear strain vs. shear stress behavior, 
Eq. (92) becomes simply 


_ sin (2 <p c ) X c 

2G\2 


(93) 


Assuming small angle approximations, a simpler expression for Eq. (93) can 
be obtained, 

= '££-■ (94) 

<-712 

Either way, the initial misalignment angle comes then as 

cp° = cp c - 7 i C - 2- (95) 

where 7 { m2 m can be defined through Eqs. (92), (93) or (94). 

Writing Eq. (92) in the form 

fcL (71-2- ) = \ sin [ 2 {<P° + 7i- 2 -)] X c (96) 


and plotting the left and right hand side of it in a (7, r) space and provides 
insight into the meaning of the solution of Eq. (92). Fig. 11(a) represents both 
sides of Eq. (96) for a material with linear shear behaviour. 


The Left Hand Side (LHS) of Eq. (96) is the shear strain vs. shear stress ma- 
terial law. The Right Hand Side (RHS) of Eq. (96) represents the shear stress 
resulting from the compressive longitudinal loading, in a rotated coordinate 
system. As the compressive stress increases, the ‘RHS’ curve corresponding 
to a general stress level shifts up, and the intersection with the ‘LHS’ curve 
defines the strain in the misalignment frame 7 lm2m . At failure (when the com- 
pressive stress equals X c ), the strain in the misalignment frame is defined as 
7i m 2 m an d is shown in Fig. 11(a). 
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(a) (b) (c) 

Fig. 11. Left and Right hand side of Eq. (96), for a material with a (a) linear 
shear behavior, (b) non-linear shear behavior, and failure by matrix cracking, and 
(c) non-linear shear behavior, and failure by instability 

For a material with non-linear shear behaviour, there could be more than 
one intersection point for each stress level, as represented in Fig. 11(b) for the 
particular case of failure onset. For such a non-linear material, the equilibrium 
position, which defines the orientation of the misaligned frame, corresponds 
to the first intersection (lower energy). As the compressive loading increases, 
the ‘RHS’ curve corresponding to a general stress level shifts up, and the 
intersection with the ‘LHS’ curve defines the strain in the misalignment frame 
7 lm 2 m. At failure (when the compressive stress equals X c ), the strain in the 
misalignment frame is defined as 7 £ m2 m and is shown in Fig. 11(b). Therefore, 
7 ±m 2 m has to correspond to the first intersection of the two curves. However, 
this is not guaranteed by the solution of Eqs. (89), (92) and (95). If the referred 
solution corresponds to the second intersection, then the solution is not valid 
and failure is due to a different mechanism, which is now discussed. 

Case 2 Kinking for pure compression as the result of instability 

A second mechanism that can promote fiber kinking is elastic instability of the 
matrix, due to the softening character of the constitutive law. As a composite 
is progressively loaded in compression, the curve ‘RHS’ in Fig. 11(b) shifts 
up, also progressively. Suppose that at the moment the two curves (‘LHS’ 
and ‘RHS’) are tangent to each other, Fig. 11(c), the matrix compression 
failure criterion (LaRC#2, Eq. (81)) is not yet verified. Then, a small increase 
in the compressive load results in the two curves not touching each other. 
Physically, this means that there is no equilibrium position and catastrophic 
failure results, due to unstable rotation of the fibers. It is concluded that 
the compressive strength measured for pure axial compression is in this case 
actually a failure due to instability, rather than matrix failure. 

The values of ip° and 7 )A 2 m corresponding to this type of failure can be ob- 
tained from the system that results from Eq. (96) and the condition that 
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expresses that the left and right hand side of Eq. (96) have the same slope at 


7l m 2" 


/ci(7l™2 m 
dfcL{llm 2 m 

<9q 


X c 


sin [2 (( p° 


= X C COS [2 (( p° + 7 lm 2 m) . 


(97) 


To summarize, the three variables p°, 7im 2 m and p c can be determined by 
Eqs. (89), (92) and (95) (matrix compressive failure) or by Eqs. (95) and (97) 
(elastic instability). 

If the instability solution occurs for Fm < 1 (in LaRC#2, Eq. (81)), then 
the instability solution must be considered; otherwise the matrix compressive 
failure solution is considered. Either way, p°, m2 m and p c are always defined. 

The initial misalignment angle tp° is a material property, and could be regarded 
as an equivalent angle that embodies microstructural defects (that can trigger 
kink-band formation) other than the initial misalignment, like oscillations in 
the fiber volume fraction or in the bonding to the resin. Knowing p° allows 
the establishment of an equation defining the misalignment angle for a generic 
plane stress situation, p : by using the transformation Eqs. (86), 

(J ii — (J 99 

fcL{hm 2 m ) = 2 — - sin (2 (p° + 7i"*2"*))+|Ti 2 | cos (2 (p° + 7 lm2m )) (98) 

and the angle p is obtained from 

F = + 7i"*2"0- (99) 

P12I 

Note that, in Eq. (98), |ti 2 | was used instead of ri 2 because it is the easiest 
way of considering simultaneously the possibility of an initial misalignment 

±p°. 

Solving this equation (where for any load case the only unknown variable is 
7 lm2 m) for non-linear shear responses is probably the most complex task of 
the whole model. For most practical cases, Eq. (98) can be simplified without 
significant error by assuming small angle approximations: 


fcL (7l m 2 m ) ~ (p° + 7l™2"0 ( — a ll + ^ 22 ) + 1^12 1 • (100) 


Solving Eq. (100) depends on the particular form of the constitutive law r = 
/cl (7)- For a linear shear behaviour, Eq. (100) was solved by Davila et al. [5], 
resulting in 


7i™ 2 ™ 


F°Gi2 + |ti 2 | 
G 12 + CTn — <7-22 


( 101 ) 
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Fig. 12. Failure envelopes for longitudinal compression acting with inplane shear for 
(a) linear shear behaviour and (b) non-linear shear behaviour 

For non-linear shear behaviour, there might be no easy way of solving Eq. 
(98) or (100) without iterating. Furthermore, Eq. (98) does not always have 
a solution, since failure by instability is also possible. If, for a specific load 
state, Eq. (98) does not have a solution (this can be easily checked by plotting 
the left and right hand side of the equation in a (r, 7) space), then failure has 
taken place by instability. The envelope for failure by instability is defined 
by the system that results from the following two conditions: (i) Eq. (98) is 
verified; and (ii) the slope (in a (r, 7) space) of the left hand side of Eq. (98) 
is the same as the slope of the right hand side. Mathematically, the following 
system follows, 

^"ll — ^"22 

sin (2 ( ip° + 7 1 m 2 m ) ) + | F1.2 1 cos (2 (9? 0 + 7 V n 2 m)) 

((Til - ^ 22 ) COS (2 (ip° + 7 i m 2 m)) - 2 |Ti 2 | Sill (2 (ip° + 7 i m 2 m)) • 

( 102 ) 

As long as no instability takes place, the misalignment frame can be defined, 
the stresses can be rotated to that frame, and a matrix failure criterion can 
be used to check for possible kink-band formation or matrix failure. 

Fig. 12 shows the application of this analysis to a biaxial compression in the 
fiber and matrix direction for a E-glass/DY063 epoxy used in the WWFE [16]. 
The material properties used are given by Soden et al. [60] and presented in 
Table 2. For the non-linear behaviour, the experimental data stress vs. strain 
points given by Soden et al. [60] are used directly by the model (Fig. 13). 

Fig. 12(a) shows the different envelopes obtained when using the LaRC04 
matrix failure criterion (Eq. (81)) and the Puck matrix failure criterion (Eq. 
(74)), considering a linear shear behaviour. 

Fig. 12(b) presents the application of this failure model for a non-linear shear 
behaviour (see Eqs. (92) and (98)). Shear non-linearity was not considered 
explicitly; instead, interpolation and extrapolation was used to get the re- 


! /ci/(7l m 2 m ) — 
dfcL (bl m 2 m ) _ 

c?7i m^m 
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Fig. 13. Linear and non-linear shear curves for an E-glass/DY063 epoxy used in the 
WWFE, experimental points from Soden et al. [60] 

quired information from the experimental data points. The linear and non- 
linear shear strain vs. shear stress curves are represented in Fig. 13. The effect 
of the non-linearity on the envelope — Fig. 12(b) — is seen to be considerable, 
which reinforces the importance of a more comprehensive characterization of 
composite materials under shear loading. 

Note the interesting implications of this model: the failure envelope is de- 
pendent on the elastic properties of the material. This is common in fiber 
compression failure models and in this case reflects the fact that failure takes 
place in a rotating misalignment frame. Indeed, the magnitude of the rotation 
of the fibers, 7 1 m 2 m, depends on the shear response. 


4-2.2 3D kinking model 

The mechanics of kink-band formation is particularly complex. Therefore, a 
physically-based model that simulates all the details of kink-band formation — 
initiation, subsequent fiber rotation, fiber fracture at kink-band boundaries, 
kink-band width, angle, broadening, and that takes into account magnitude, 
distribution and orientation of initial defects, in-plane shear non-linearity, and 
all this for a generic 3D stress state — does not seem possible. For the validation 
of emerging, simplified models, an added difficulty results from the difficulties 
in obtaining trustworthy experimental data for compression failure under a 
multiaxial stress state. Consequently, the validation of a specific mechanical 
model for fiber kinking is bounded to some controversy. 

The 2D model for fiber kinking just presented assumes that the initiation of 
kink-band formation is triggered by matrix failure due to initial fiber misalign- 
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Fig. 14. 3D kinking model 

ment. As with most fiber-kinking models, the 2D model assumes that kinking 
happens in the plane of the lamina. 

However, there is significant evidence supporting the importance of 3D anal- 
yses for understanding failure by kink-band formation [61,67-71,73-78]. The 
model that follows generalizes the previous one for a generic 3D stress state. 

Consider a unidirectional lamina under compressive stress as shown in Fig. 
14(a). The stresses acting on the (2,3) plane are shown in Fig. 14(b). The 
kink plane is at an angle ip with the 2 axis, as shown in Figs. 14(c) and (d). 
The actual value of the angle ip depends on the particular stress state 3 . A 2D 
kinking model assumes that the angle ip in Fig. 14(c) and (d) is equal to zero. 
If the composite is constrained so that it cannot move laterally, then the kink 
plane would have an angle ip = 90°. In general, ip will have a value between 0 
and 180°. 

The local stresses in the (2, 3) plane can be found as a function of ip through 
the use of transformation equations. For a potential kinking plane oriented 
at an angle ip , a negative stress a 3 ^ will act as to close microcracks in the 

3 In reality, defects as fiber initial misalignment are not necessarily homogeneously 
distributed and the kinking plane orientation could in fact also be influenced by a 
particular distribution 
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matrix for alternative kinking planes, while a positive cr 2 ^ 2 ^ will act as to 
open microcracks thus favouring kinking in that plane. Therefore, the kinking 
plane is expected to be created at the orientation that maximizes cr 2 ^ 2 ^ and 
minimizes which coincide with the local principal directions. Another 

argument is that the reduced shear stiffness in the kink band would result 
in fibers rotating in different directions, in case r 2 v> 3^ was nonzero — however, 
experimental observations support the contention that kink bands lie in a 
plane. 

With the assumption that the kink plane happens at an angle such that 
r 2 i, 3^ = 0, the value of the angle if} that defines the kink plane is given by 


tan (2-0) 


2^23 

&22 ~ 0~33 


(103) 


The stresses in the kink plane are defined as 


a 22 + &33 <7 22 — C33 / n , X | . / 0 , \ 

c x 2 V-2V- = t cos (2 Ip) + T -23 sin (2ip) 

&3'l>3 = cr 22 + °"33 — Cr 2 i/>2’/> 

< r 12i , = T 12 cos (^) + T31 sin (-if)) 

T 2 i/>3V> = 0 

T3V-1 = T31 cos (^) - n 2 sin (^) . 


(104) 


After defining the kink plane, the stresses are then rotated to the misalign- 
ment frame. This misalignment frame is defined by first determining 7 lm2?n by 
solving the following equation iteratively: 


fcL (7i m 2 m ) — 2 s ^ n d - 7i m 2 m )) T |ri 2 v,| cos (2 (93 + 7l m 2 m ) ) • 

(105) 

After obtaining 7 lm2 m from Eq. (105), the angle <p is obtained from 

A = ((p° + 7lm 2 m) • (106) 

|Tl2*| 

If Eq. (105) has no solution, then failure has taken place by instability. Oth- 
erwise, the misalignment frame is defined and the stresses can be rotated to 
it to check for matrix failure. In the misalignment frame, the stresses can be 
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written as 


(7 ^ 771^71 


(Tli + <7 2 020 , 

2 + 




cos [Zip) + r 12 0 


( 7 2 rn 2 rn — 0*11 “f" CT 2 7 /> 2 ^ — ( 7 \m-^m 

( 7 11 — (7 2 ^ 2 ^ . / 0 X / 0 \ 

< rim 2 m = sin (2<p) + r 12 </> cos 

t 27773^ — 'T’2'03'0 cos ( 92) 7-37/, 1 sin (9^) 

7 ^ip-^m — — T3 COS (9^) • 


( 107 ) 


The LaRC04 criteria for matrix tension failure, Eq. (40), and matrix compres- 
sion failure, Eq. (81), can be applied in the misalignment frame to predict 
failure. The type of failure predicted is matrix failure, and this might or not 
promote subsequent fiber kinking. While it might be trivial in some cases to 
predict that subsequent fiber kinking takes place — like for pure compression 
in the fiber direction, the same is not true for all load combinations. 


For <T 2 m 2 m < 0, it seems reasonable to assume that failure is by kink-band 
formation only if Eq. (81) is verified with a = 0. Otherwise, the failure is 
considered to lead to matrix failure, without kink-band formation. Thus, kink- 
band formation with cr 2 m 2 m < 0 is predicted with the criterion 

LaRC04 #4 FI f = oT < 1 (108) 

bis - V L (T 2 m 2 m 

while matrix failure under biaxial compression is predicted by 

/ Tm \ 2 / Lm \ 2 

LaRC04 #5 FI M = y S T_ rj T a m J + \s L - n L v™ ) ~ 1 (1 ° 9) 


with 


m ( J 2 m 2 m + <73031/1 0 2 m 2 m — <T 3 0 3 0 ■ / r > \ 

ct" = 1 cos { 2 a) + r 2 m 3 0 sin { 2 a) 


n 


r Tm _ 


& 2 ™ 2 ™ ~ <73031/1 ■ / n \ , (n \ 

sm {2a) + r 2 m 30 cos {2a) 


( 110 ) 


r Lm = Tim 2 m cos (ex) + 7303m sin (a) 


where the angle a, which is comprised in the interval ]0, 7r[, is obtained by 
trying a small number of tentative angles. 


For cr 2 m 2 m > 0, it seems more difficult to agree on a criterion for eventual fiber- 
kinking after matrix failure, in the absence of experimental data. Possible 
solutions to identify the conditions triggering fibre-kinking are the use of a 
threshold value for <7n, or to consider fibre-kinking when the last term of Eq. 
(40) is equal to 1. Without further support from experimental data, X c /2 is 
taken here as the threshold value. Thus, for <r 2 m 2 m > 0 and from Eq. (40), 
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the criterion for matrix tensile failure under longitudinal compression (with 
eventual fiber-kinking) is, 


LaRC04 #6 


fIm/f — (i — g) 
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5 Applications 


5.1 Failure envelope (a 22 , t 12 ) 


A comparison of results from various failure criteria with the experimen- 
tal results in the (<t 2 2 , t 12 ) stress plane for an unidirectional composite E- 
glass/LY556, is shown in Fig. 15. The material properties and test results 
were reported in the WWFE by Soden et al. [60,86]. The compressive, tensile 
and shear strengths shown in Table 3 were obtained by averaging the corre- 
sponding test data presented in Refs. [60,86]. To apply the Puck tensile matrix 
failure criterion, it was assumed that (dri 2 / da 2 2 )\ a2=0 - = {dr\ 2 /da 22 ) | 0+ . 

The envelopes for Hashin’s 1980 criteria and Sun’s criteria were calculated us- 
ing a transverse shear strength obtained from Eq. (68). The friction parameter 
r] L for Sun’s criteria was obtained from Eq. (75) with r] T from Eq. (67). 

It can be observed that within the positive range of cr 22 , all the quadratic 
failure criteria and LaRC04 give satisfactory results. The most interesting 
behavior develops when cr 22 becomes compressive. Hashin’s 1973 criterion gives 
an elliptical envelope with diminishing ti 2 as the absolute value of compressive 
022 increases, while the experimental data shows a definite trend of shear 
strength increase as 022 goes into compression. 

The envelope for Hashin’s 1980 criteria provides a modest improvement in 
accuracy compared to the 1973 criterion. Of the criteria shown in Fig. 15, 
Sun (Eq. (76)), LaRC04 #2 (Eq. (81)), and Puck (Eq. (74)) capture the shear 
strength increase at the initial stage of compressive 022 • The results indicate 
a significant improvement over Hashin’s criteria. 

Fig. 16 shows the application of LaRC04 for AS4/55A, tested by Swanson et 
al. [87]. The material properties used are given in Table 3. The predictions 
for T800/3900-2 composite tested by Swanson et al. [77] are shown in Fig. 17, 
and the material properties in Table 3. 
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Table 3 

Mechanical properties of E-glass/MY750, from Soden et al. [60,86] 


Material 

En 

E 22 

G 12 

S L 

Y t 

Y c 

Ot 0 


(GPa) 

(GPa) 

(GPa) 

(MPa) 

(MPa) 

(MPa) 

(°) 

E-Glass/LY556 

53.48 

17.7 

5.83 

66.5 

37.5 

130.3 

53 

AS4/55A 

126 

11 

6.6 

51.3 

27 

91.8 

53 

T800/3900-2 

155 

8.5 

5.5 

100.9 

48.8 

201.7 

53 


x 
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Fig. 15. Failure envelopes and WWFE test data for unidirectional composite 
E-Glass/LY556 

5.2 Failure envelope (an, < 722 ) 



All of the failure modes represented by the six LaRC04 criteria (Eqs. (40), 
(81), (82), (108), (109) and (111), and summarized in the Appendix) can 
be represented as a failure envelope in the (<J \ , a 2 ) plane. Testing for biaxial 
loads presents a number of complexities, and experimental results are rare. 
However, Waas and Schultheisz [88] report a number of references in which 
multiaxial compression was studied by superposing a hydrostatic pressure in 
addition to the compressive axial loading. For all materials considered, there 
is a significant increase in compressive axial strength with increasing pressure. 

As an example, a test case from the WWFE [1] is studied here: the biaxial 
compression of a 0° E-glass/MY750 epoxy lamina, with material properties 
given in Table 4. In addition to those properties, the in-plane Poisson’s ratio 
is z/ 12 = 0.278 and the fracture angle for pure in-plane transverse compression 
is a D = 53°. For the application of the model with non-linear shear behaviour, 
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Fig. 16. Failure envelopes and test data from Swanson et al. [87], for an AS4/55A 
composite 



O 22 (MPa) 

Fig. 17. Failure envelopes and test data from Swanson et al. [77], for an T800/3900-2 
composite 

the experimental curve given by Soden et al. [60] is used directly. The linear 
and non-linear shear laws are plotted in Fig. 18(a). 

For the LaRC04 3D kinking model, the kinking angle is predicted using Eq. 
(103). The magnitude of <73 is unknown, and depends on the existence of 
constrains on the 3 direction. Therefore, the 2 roots of Eq. (103) must be 
considered and are ipx = 0° and = 90°, meaning that the kink plane can 
either be in the plane of the lamina, or in the through-the-thickness direction. 
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Table 4 

Mechanical properties 


Longitudinal Transverse Shear 

Tensile Compression Tensile Compression 

Modulus (MPa) 1280 800 40 145 73 



(a) (b) 

Fig. 18. (a) Linear and non-linear shear law; (b) biaxial (a \ . < 72 ) failure envelope of 
0° E-glass/MY750 epoxy lamina, assuming a kink band in the plane of the lamina 

Assuming first that the kink band develops in the plane of the lamina (0 = 0°), 
either due to the micromechanics of the material or imposed by the testing, the 
failure envelope comes as in Fig. 18(b). Puck’s analysis results, which showed 
a particularly good correlation with experimental results in the WWFE [1], 
is also shown for comparison. In the figure, there is good agreement between 
LaRC04 and Puck in all quadrants except biaxial compression, where LaRC04 
predicts — for a kink band developing in the plane of the lamina — an increase 
of the axial compressive strength with increasing transverse compression. 

If it is assumed that the kink band is formed in the through-the-thickness 
direction (second root of Eq. (103), 0 = 90°), the envelope shown in Fig. 
19(a) is predicted. If the orientation of the kink plane is unrestricted, the 
envelope in Fig. 19(b) is obtained. 


5.3 Axial compression with superposed hydrostatic pressure 


There is not much experimental data on fiber kinking under a multi-axial 
stress state. One exception is the compressive behaviour of composite rods 
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Oil 

(MPa) 



Fig. 19. Biaxial (< 711 , 022 ) failure envelope of 0° E-glass/MY750 epoxy lamina (a) 
assuming through-the-thickness kinking; (b) assuming that there are no restrictions 
to the kinking plane 

with superposed hydrostatic pressure [89,90]. 

Wronsky and Parry [89] made a characterization of the effect of hydrostatic 
pressure on the compressive strength for a glass fiber reinforced composite. The 
compressive strength without superposed hydrostatic pressure was measured 
as 1150 MPa. Three values of shear strength are reported, depending on the 
test method and specimen dimensions, 42, 48 and 59 MPa. The biggest value of 
the three is used herein, because it is expected to be the most representative of 
the three. These material properties are not sufficient to completely define the 
material. Therefore, some material properties have to be assumed from their 
typical values. The shear modulus is 6.6 GPa, the fracture angle in matrix 
compression is a Q = 53°, and the transverse compressive strength is 140 MPa. 
Linear and a non-linear approximations of the material’s shear response are 
shown in Fig. 20(a). The non-linear curve uses the logarithmic 4 function r = 
k\ In (/c 2 T T 1 ), with ki = 200 MPa and k 2 = 33, which yields the correct 
initial shear modulus, G\ 2 = k\k 2 = 6.6 GPa, and matches particularly well 
the experimental points, as shown in Fig. 20(b). The same figure shows that, 
considering a linear shear model, the global trend is captured, but an over- 
prediction of the axial strength results for the higher values of superposed 
hydrostatic pressure. 

Parry and Wronsky [90] compressed 8 mm carbon-epoxy rods for several values 
of superposed hydrostatic pressure. Three values of compressive strength at 
atmospheric pressure are reported (depending on the test specimen’s batch 
and its shape): 1.3 MPa, 1.33 MPa and 1.5 MPa. The corresponding average 
is considered here. The shear strength reported is 75 MPa. As before, some 
material properties have to be assumed to apply the present model. The shear 

4 A logarithmic function is used because, while correctly representing the shear 
non-linear behaviour of composites, it can be easily inverted 


43 



Fig. 20. Test case corresponding to Wronsky and Parry’s [89] experiments (a) linear 
and non-linear shear curves considered; (b) compressive strength as a function of 
the hydrostatic pressure 




(a) 


(b) 


Fig. 21. Test case corresponding to Parry and Wronsky ’s [90] experiments (a) linear 
and non-linear shear curves considered; (b) compressive strength as a function of 
the hydrostatic pressure 

modulus is 6.6 GPa, the fracture angle in matrix compression is a Q = 53°, the 
transverse compressive strength is 170 MPa. The predictions for a linear shear 
model (Fig. 21(a)) are shown in Fig. 21(b). The effect of shear non-linearity 
is the same as for the previous example. Considering the non-linear stress 
vs. strain curve in Fig. 21(a) (logarithmic curve, r = k± In (fc 2 7 + 1), with 
ki = 120 MPa and = 55, which yields the same initial shear modulus), the 
predictions in Fig. 21(b) are in good agreement with the experimental points. 

The comparison of experimental results and predictions in Figs. 20 and 21 
suggests that the physics of the compressive behaviour is captured by the 
kinking model. 


5-4 Cross-ply laminates 


Sliuart [73] studied the compression failure of [±0] s laminates and found that 
for 6 < 15°, the dominant failure mode in these laminates is interlaminar 
shearing; for 15° < 6 < 50°, it is in-plane matrix shearing; and for 9 > 50°, it 
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Table 5 

Material properties of AS4/3502 


An 

E 22 

G 12 

"12 Sfi 

Y c 

A c 

(GPa) 

(GPa) 

(GPa) 

(MPa) 

(MPa) 

(MPa) 

127.6 

11.3 

6.0 

0.3 95.1 

255 

1045 


Compressive strength (MPa) 



Lamination angle (°) 


Fig. 22. Compressive strength as a function of ply orientation for AS4-3502 [±0] s 
laminates 

is matrix compression. Fiber scissoring due to matrix material non-linearity 
caused the switch in failure mode from in-plane matrix shearing to matrix 
compression failure at larger lamination angles. The material properties used 
for the analysis are given by Shuart [73] and shown in Table 5. Further to those 
properties, the fracture angle in pure transverse compression is considered to 
be a Q = 53°. 

The results shown in Fig. 22 indicate that the predicted strengths using 
LaRC04 criteria correlate well with the experimental results. The compres- 
sive strength predicted using Hashin’s 1980 criteria is also shown in Fig. 22. 
For 9 < 20°, the Hasliin criterion for fiber compression results in an over- 
prediction of the failure load because the criterion does not account for the 
effect of in-plane shear on fiber failure. For lamination angles near 70°, the 
use of the Hashin criteria result in an under-prediction of the failure load 
because the criteria do not account for the increase in shear strength caused 
by transverse compression. All of these effects are correctly captured by the 
LaRC04 criteria, which results in a good correlation between the calculated 
and experimental values. 
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6 Conclusions 


Material models based on the physics of failure can yield criteria that succeed 
in predicting failure under a broad range of load combinations. Consider- 
ing the current difficulties in accurately predicting failure [1], the ability of 
physically-based models to capture the physics of the failure process and cor- 
rectly predict failure envelopes is particularly significant and should encourage 
further research. In this paper, physically-based models of the different failure 
modes result in a set of six physically-based failure criteria devoid of empirical 
variables — the LaRC04 criteria — that can be used directly by the designer. A 
few considerations on the models presented are now discussed; some of them 
will be the focus of further work. 

Non-linear shear behaviour is an important factor affecting strength in a gen- 
eral loading situation. Furthermore, it is known that the shear behaviour de- 
pends on the stress in the in-plane transverse direction (the shear response 
gets “stiffer” with compressive transverse stresses [6]). In the future, it would 
be desirable to consider the dependence of the shear stress on eventual su- 
perposed transverse stress. Furthermore, the shear strength reported in the 
literature is often the shear stress at some strain level — as recommended by 
actual standard methods [91]. Although convenient for design purposes, this 
value of the shear stress has no direct correspondence with failure and is mean- 
ingless for physically-based failure models. 

For fiber kinking, only matrix failure (as well as some form of instability) is 
assumed to trigger the failure. This consideration is believed to be accurate for 
carbon-reinforced composites. However, a more general model that includes 
the possibility of fiber buckling, as well as fiber failure itself would be desirable. 
Furthermore, the strength in a generic loading situation is very sensitive to the 
initial fiber misalignment angle. However, this misalignment angle results from 
manufacturing limitations, and should be expected to be different for panels 
or components with different dimensions and geometry. Therefore, the value 
measured using standard compression specimens could not be representative 
for other more complex components. Finally, residual stresses — present in the 
compression specimens — affect the value measured of the fiber misalignment 
angle. In more complex components, the residual stresses will certainly be 
different, and relax over time, which restricts the applicability of the proposed 
model, and should be addressed in the future. 
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Appendix - Pragmatic solution for the critical energy release rates 
for non-linear shear behaviour 


Introduction 


Consider a generic solid made of a linear elastic material with an initial crack — 
Fig. 23(a). Let the crack grow by an area d A. The energy absorbed by the 
solid to create the surface of area dA is dW s . Neglecting thermal effects, energy 
balance requires d W s to be equal to the energy transferred to the body through 
external work, dF, minus the change in strain energy dU of the body, i.e., 


dW s = d F - dU. 

By definition of energy release rate G, it follows that 

dW s _ dF — dU _ dU 
T “ ~dA ~ dA “ ~dA 

where II is the potential energy defined as II = U — F. 


( 112 ) 


(113) 


Considering the diagram in Fig. 23(b), it can be concluded that 


dU = ^ (P + dP) (A + dA) - ^PA 


and 


(114) 



Fig. 23. (a) Generic cracked body being loaded; (b) load displacement curve corre- 
sponding to a crack length 2a (i) and 2 (a + da) (ii) 
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(115) 


Thus, 


dF = Pd A + -dPdA 
2 



8A 


d(F-U) = ^ (P<9A - A<9P) 
and the critical energy release rate is expressed as 

r = 1 A— - API A— _ Al 

' 2 f 9,4 <V. l J 25 V da da ) 


(116) 


(117) 


where b is the thickness of the solid and a is the crack length. Defining the 
compliance C of the cracked body as 


C = 


A 

P 


(118) 


it follows that 

dA = PdC + CdP (119) 


and, designating the load at onset of crack propagation as 
energy release rate is 


G c 


P c 2 dC 
2b da ' 


Pr 


the critical 

( 120 ) 


Non-linear in-plane shear behaviour 


Considering the energy release rate given by Eq. (113), but carrying the differ- 
entiation at constant displacement of the load application-points, hence with 
no work done by the external forces, it follows, at onset of crack propagation, 



= G c . 


( 121 ) 


Following the approach first proposed by Rivlin and Thomas [92] , the assump- 
tion is made that the presence of a crack in a body will reduce the strain-energy 
density to zero over a well-defined volume V, and will not affect it outside that 
area. Let the uniform strain density of the uncracked body be 

W Q = f (JijdEij ( 122 ) 

J 0 

and for the particular case of pure in-plane shear loading 

T712 

W 0 = / r 12 d7i2. (123) 

Jo 
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The loss in strain energy due to the presence of the crack is 


-{U- U 0 ) = W 0 V = (712) v (124) 

where 

/■712 

X (712) = 2 / ^12^712 (125) 

Jo 

and U Q is the strain energy of the body for the situation without a crack. The 
critical energy release rate is 

G//c=)x(7iy^. (126) 

Consider first the situation where the elliptical crack grows in the transversal 
direction. The area of the crack is proportional to cl 0 cil and thus 

d A oc a L da 0 . (127) 


The volume where the strain energy density is reduced to zero due to the crack 
can be assumed to be proportional to a 2 Q aL and thus 


dV oc a 0 aLda 0 . 


It follows from Eqs. 127 and 128 that 


dV 

dA 


oc a 0 - 


(128) 


(129) 


In order to recast Dvorak and Laws’s solution for linear clastic materials [9], 
the constant of proportionality in Eq. 129 has to be n. Thus, the expression 
for GjiciT ) results: 

7TC1 

Gnc(T ) = -f X (7^2) • (130) 

Proceeding in the same way for propagation in the longitudinal direction, the 
following expression results: 


G IIc (L) 


it a 0 


X (in) ■ 


(131) 
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Appendix - Summary of LaRC04 failure criteria 


The LaRC04 criteria is summarized in Table 6. The required unidirectional 
material properties for the criteria are: E u , E 2 2 , G 12, v- 12 , X T , X c , Y' r , Y c , 
S L , Gic, Gji c . The following optional properties may be provided: i] L , a Q - The 
expressions needed to obtain the parameters used in in Table 6 are given in 
following sub-sections. 


Expressions for matrix tensile failure 


The in-situ strength values for thick embedded plies are 

Y? = 1.12 V2 Y t and 

Y ixis = X " 1 [ 2 x (Y12)] • 

The in-situ strengths for thin embedded plies are: 


(132) 


Y t = 

IS 


8 G 


Ic 


TCtA° 


22 


where A 22 = 2 


Ti2iis = X 1 


8G 


lie 


Tit 


'21 


E- 


22 


E 


11 , 


(133) 


The toughness ratio g is obtained from fracture mechanics test data or from 
the unidirectional properties. For a thin embedded ply, 


9 = 


G 


Ic 


G 


(134) 


lie 


otherwise, g can be obtain from either of the following expressions: 

■Y-Y ) 2 


g = 1-12 


x(i 


(135) 


12 ) 
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Expressions for matrix compressive failure 


The fracture angle for pure transverse compression can be considered as a 0 = 
53° if no experimental value is available. The friction parameter r] T is obtained 
from 


tan(2a 0 ) = — 
V 


(136) 


the transverse strength S T is defined as 


S T 


Y° cos (a 0 ) I sin (a G ) + 


cos (a 0 ) 
tan (2a 0 ) 


(137) 


and, in the absence of experimental data to obtain r) L , this is defined from the 
relation 

T T 

T] T) 

sh = s 


(138) 


The stresses in the potential fracture planes (in the interval a G [0, 7r[) are 
obtained with 


<722 + 033 . 022-033 / n \ , ■ f n \ 

0 n. = z t cos (2a) + r 23 sin [ 2 a ) 


< T 1 = 


$■ 


022 + 033 022 — 033 


sin (2a) + r 2 3 cos (2a) 


(139) 


t l = r 12 cos (a) + r 3 i sin (a) . 


Expressions for fiber kinking and for biaxial compressive failure 


For the determination of 7^ m2 m, <p° and (p c , a two steps solution is required, 
since kinking can be promoted by instability or matrix failure 5 . The solution 
for instability results from solving the following system for ip° and 7J m2m 


/cl( 7r»: 

dfcU 7i m 2171 ) 

im 2171 


'Jin 


x c 

--^- Sill [ 2 (<P° + 

= —X c cos [2 + 7i m 2 m ) ] 


(140) 


and ip c follows as <p c = <p° + 7^ m 2 m. If this solution is verified with F M < 1 
(in LaRC04 Eq. (81) or from Table 6), then the instability solution must 
be considered; otherwise the matrix compressive failure solution is considered, 
and <p c , 7i m 2 m and <p° are defined as follows: 


5 For a linear material, there is no need to check for the instability solution 
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ip c = arctan 


1- \ 

1 

T— 1 

-4| 

(ft + v L ) 

b 

21 

(ft +ft) 

1 


7i»2- = fcl () sil > ( 2 7> c ) V c ) 


(141) 

(142) 


where /cl is the in-plane shear constitutive shear law. For a linear behaviour 
it simplifies to 

ip c X c 


7l™2” 1 


G 


12 


Either way, the initial misalignment angle comes then as 


C P — • 


The plane where kinking takes place is defined by 

tan (2ip) = — — — — 

@22 ^33 

and the stresses rotated to that plane are 

@22 + @33 , @22 ~ @33 


(143) 

(144) 


(145) 


@2^2'!’ — 


+ 


2 2 

@yl>yl> = @22 @33 @2^2^ 

r 12 i> = T 12 cos (^) + t 3 1 sin (^) 

@ 2 ^ 3 ^ = 0 

T3V-1 = T31 cos (ip) - T12 sin (^) . 


cos (2ip) + r 2 3 sill (2ip) 


(146) 


The shear strain in the misalignment frame 7 “ m2 ™. is obtained by iteratively 
solving the equation 


/ci( 7 i™ 2 m ) ~ 2 sin (2 (</? + 7 i m 2 m )) T l' r i 2 V’| cos (2 ( 9 ? + 7 i m 2 m )) 

(IV) 

which, for a linear shear behaviour and with small angle approximation reduces 
to 


7l m 2 m 

and the angle tp comes as 


p >° G \2 + |t 12 | 
G 12 + @11 ~ a 2 2 


(148) 


<P= I |(^° + 7l-2-)- ( 149 ) 

Fl2*|. 

Eq. (147) does not always have a solution, since failure by instability is also 
possible. If, for a specific load state, Eq. (147) does not have a solution (this 
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can be easily checked by plotting the left and right hand side of the equation 
in a (7, r) space), then failure has taken place by instability. The envelope for 
failure by instability is defined by the following system: 

/ci/(7i m 2 m ) = 2 sin (2 {ip + 7i m 2 m )) ~k P12 1 /’ I cos (2 yp + 7i m 2 m )) 

dfcUi^'u) = _ ( aii _ ^ cos ("2 (93 0 + 7 lm2m )) : - 2 |t 12 ^| sin (2 (ip° + 7 1 m 2 m)) . 

(150) 

If Eq. (147) has a solution (and thus the material has not failed by insta- 
bility), matrix failure is checked next. After knowing the orientation of the 
misalignment frame, the stresses can be rotated to it using 


OH + (72^2^ . On — O 2^2'!’ \ , . / 0 \ 

Oim lm = 1 cos [2(f) + t 12 ip sm [2<p) 


U2m2m — Oil T O 2i/’2V’ — Olm\ 

Oil — o 2^2^ ■ /r 

T |r«2 m — sin 


sin (2 tp) + t 12 v< cos (2 p) 


T 2 m 3 4 , — T2i/ , 3’/’ cos {p) 'Ti’/’i sin {p) 

T^\m = T3V'iV' COS (92) • 


For biaxial compression failure, the rotation to the potential fracture planes 
(in the interval a G ]0, 7r[) are obtained with 



(T2">2 m + 0‘3’/'3’/' . 0‘2m 2 m — (J3V.3V- / 0 \ . N 

1 — — cos {2a) + r 2 m 3 v; sm {2a) 


O 2™2™ — 0 3 rp 3 ^ . , n \ , tn \ 

sm {2a) + T 2 m 3 ^ cos {2a) 

Tim 2 m cos (cc) + r 3 v>im sin (a) . 


(152) 


For <7 2 m 2 m > 0, it is considered that hber kinking follows matrix failure only 
if a n < —X c /2. 
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Appendix - Particularization of the non-linear shear behaviour for 
a 3 r</ order polynomial 

The LaRC04 criteria presented here accounts for a generic non-linear shear be- 
haviour. Any expression for non-linear shear can be used, such as a polynomial, 
a logarithmic, any other best fit to experimental data or the experimental data 
itself. This appendix presents the particularization of the criteria for Hahn and 
Tsai’s 3 rd order polynomial representation for shear non-linearity [50]: 

712 = 77 -T 12 + (3 ( ti 2 ) 3 (153) 

Lri2 

with (3 obtained from the shear tests (/? = 0 for linear behaviour). 
Expressions for matrix tensile failure 


The in-plane shear strength of thick embedded plies and thin plies have the 
same general form, 



where the parameter <f> is obtained as: 



for a thick embedded ply (155) 

for a thin ply (156) 


A linear shear behaviour is obtained when the constant (3 tends to zero in Eq. 
(153) 



The particular values of the in-situ shear strengths for a linear shear behaviour 
can thus be expressed as 


St = V2 S L 



for a thick embedded ply 
for a thin ply 


(158) 

(159) 
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Expressions for fiber kinking and for biaxial compression failure 


Initial misalignment angle For the determination of 7 f m2m , (p° and ip c , a 
two steps solution is required, since kinking can be promoted by instability 
or matrix failure 6 . The solution for instability results from solving iteratively 
the following system for <7° and rf m2 ™ 




r- 




Sill 


2 iy° + cL T + P ( r l-2-) 3 


G 12 A ^ a ( r l m 2 m ) 

and then compute 7)m 2 m and Lp c as 


X C COS 2 (y)° + p^-jfm 2 rn + / 3 (t^™) 3 ) 


(160) 


7l m 2 m < ^"l m 2 m “1“ fi {r yrn^rti ) 

bri2 

^ I • 


(161) 

(162) 


If this solution is verified with Fm < 1 (in LaRC04 Eq. (81) or from Table 
6), then the instability solution must be considered; otherwise the matrix 
compressive failure solution is considered, and ip c is determined from Eq. (141), 
7 {m 2 rn comes from 

7i"» 2 m = ^ Q sin ( 2( P C ) A " C ) + P Q sin ( 2< d c ) A " C ) (163) 
and (p° = <p c — 7 \m 2 m. 


Current misalignment frame The shear strain in the misalignment frame 
7 lm2 m is obtained by iteratively solving the equation 

1 i 

7i m 2 m = 71 T (di m 2 m ) T P It (7i m 2 m )] (164) 

Vti2 

with 

T*( 7 lm2 m) = 1 - 20 sin (2 (tp° + 7im 2 m)) + |t 12 V-| COS (2 (<7° + 7im 2 m)) • 

(165) 

Eq. (164) does not always have a solution, since failure by instability is also 
possible. If, for a specific load state, Eq. (164) does not have a solution (this can 
be easily checked by plotting the left and right hand side of the equation), then 

6 For a linear material, there is no need to check for the instability solution 
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failure has taken place by instability. The envelope for failure by instability is 
defined by the following system: 


T 


CTii — O r 2 1 /’2 




s i n p + + (3 (r lm2 m) )) T 


+ |ri2^| cos (^2 (</?° + ^r lm2m + a (r lm2m ) 3 )) 

1 —^7 p = - (o-n - cos (2 (V + ^r lm2m + /? (r lm2 m) 3 )) 


Gl2 


-2 |ri 2 7/,| sin (2 (</?° + ^ r lm2m + [3 (T lr a 2 m)' 


( 166 ) 


64 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


i 

The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Department of Defense, Washington Headquarters Services, Directorate for Information Operations and 
Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person 
shall be subject to any penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. 

PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 


1. REPORT DATE (DD-MM-YYYY) 2. REPORT TYPE 3. DATES COVERED (From - To) 

01- 02 -2005 Technical Memorandum 

4. TITLE AND SUBTITLE 

Failure Models and Criteria for FRP Under In-Plane or Three-Dimensional 
Stress States Including Shear Non-Linearity 

5a. CONTRACT NUMBER 

5b. GRANT NUMBER 

5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 

Pinho, Silvestre T.; Davila, Carlos G.; Camanho, Pedro P.; 
Iannucci, Lorenzo; and Robinson, Paul 

5d. PROJECT NUMBER 

5e. TASK NUMBER 

5f. WORK UNIT NUMBER 

23-064-30-22 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Langley Research Center 
Hampton, VA 23681-2199 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

L-19089 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

10. SPONSOR/MONITOR'S ACRONYM(S) 

NASA 

11. SPONSOR/MONITOR'S REPORT 
NUMBER(S) 

NASA/TM-2005-2 13530 


12. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 39 

Availability: NASA CASI (301) 621-0390 

13. SUPPLEMENTARY NOTES 

Pinho, Iannucci, and Robinson: Imperial College; Davila: Langley Research Center; Camanho: University of Porto 
An electronic version can be found at http://ntrs.nasa.gov 


14. ABSTRACT 

A set of three-dimensional failure criteria for laminated fiber-reinforced composites, denoted LaRC04, is proposed. The 
criteria are based on physical models for each failure mode and take into consideration non-linear matrix shear behaviour. The 
model for matrix compressive failure is based on the Mohr-Coulomb criterion and it predicts the fracture angle. Fiber kinking 
is triggered by an initial fiber misalignment angle and by the rotation of the fibers during compressive loading. The plane of 
fiber kinking is predicted by the model. LaRC04 consists of 6 expressions that can be used directly for design purposes. 
Several applications involving a broad range of load combinations are presented and compared to experimental data and other 
existing criteria. Predictions using LaRC04 correlate well with the experimental data, arguably better than most existing 
criteria. The good correlation seems to be attributable to the physical soundness of the underlying failure models. 

15. SUBJECT TERMS 

Failure criteria; Matrix cracks; Fiber kinking; Damage; Composite materials 


16. SECURITY CLASSIFICATION OF: 


a. REPORT b. ABSTRACT c. THIS PAGE 


17. LIMITATION OF 
ABSTRACT 


18. NUMBER 
OF 

PAGES 


19a. NAME OF RESPONSIBLE PERSON 

STI Help Desk (email: help@sti.nasa.gov) 

19b. TELEPHONE NUMBER (Include area code) 


(301) 621-0390 


Standard Form 298 (Rev. 8-98) 

Prescribed by ANSI Std. Z39.18 


U 


U 


u 


uu 


69 


































